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SUMMARY:  This  report  presents  explicit  analytical  expressions  for  the  thermal  expansion  of  a 
rod  with  heat  input  at  one  end  and  forced  convection  cooling  along  the  length  of  the  rod.  When 
the  rod  is  graduated  for  length  measurement  (i.e.,  a ruler  or  ball  screw),  thermal  expansion  causes 
an  error  in  the  length  measurement.  For  the  problem  described  here,  the  thermal  expansion  or 
length  error  depends  upon  both  distance  along  the  rod  and  time.  Calculations  are  presented 
illustrating  the  magnitude  of  the  measurement  error  for  representative  conditions  of  heating  and 
cooling. 

The  report  considers  two  heat  input  models:  step  temperature  input  and  step  heat  flux  input.  A 
representative  example  illustrates  the  system  response  to  the  two  inputs.  The  detail  calculations 
illustrate  the  differences  between  the  two  idealized  input  models.  These  idealized  models  are 
basis  response  functions  that  may  be  used  to  estimate  system  response  to  general  time-varying 
heat  input  conditions.  The  implications  for  using  these  results  for  system  metrology  are 
discussed. 
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STATEMENT  OF  THE  PROBLEM 


Ball  screws  are  a common  mechanical  system  used  to  position  an  object.  A ball  nut  is  attached 
to  the  ball  screw  and  as  the  ball  screw  is  rotated  the  ball  nut  moves  along  the  screw.  The 
distance  that  the  ball  nut  moves  with  each  rotation  is  established  by  the  constant  pitch  of  the  ball 
screw.  The  positioning  accuracy  of  this  system  is,  ideally,  dependent  upon  the  manufacturing 
precision  of  both  the  ball  screw  pitch  and  the  matching  ball  nut.  The  ball  screw  and  ball  nut  are 
manufactured  and  calibrated  to  establish  the  positioning  accuracy  of  the  system.  This  calibration 
is  conducted  under  controlled  conditions  of  temperature  and  low  speed  rotation  of  the  ball  screw 
to  ensure  calibration  accuracy.  The  above  considerations  also  apply  to  a graduated  scale  or  ruler 
although  in  this  instance  no  motion  is  involved. 

When  the  ball  screw-ball  nut  system  is  used  in  a mechanical  system  such  as  a machine  tool,  the 
operating  conditions  are  quite  different  from  the  calibrated  conditions.  Although  the  machine 
tool  may  be  independently  calibrated  to  determine  the  positioning  accuracy  of  the  system,  this 
calibration  may  not  completely  encompass  the  entire  range  of  operating  conditions  of  the 
machine.  Temperature  changes  from  the  calibration  conditions  are  commonly  recognized  as  one 
of  the  more  important  factors  affecting  positioning  accuracy  during  operation.  The  effect  of 
temperature  changes  on  the  positioning  accuracy  of  the  ball  screw-ball  nut  system  is  only  one 
aspect  of  the  more  comprehensive  (and  complicated)  problem  for  the  entire  machine. 

Considerable  design  effort  and  expense  is  expended  to  provide  a length  measurement  system  for 
a machine  tool  that  is  independent  of  the  ball  screw-ball  nut  system.  These  systems  require 
computer  control  of  the  machine  positioning  using  the  independent  length  measurement  to 
control  the  motion.  A scale-based  metrology  system  may  represent  five  percent  of  the  machine 
cost  and  a laser-based  metrology  system  may  represent  25-35  percent  of  the  machine  cost.  These 
costs  are  in  addition  to  the  computer  control  system.  Cost  considerations  commonly  result  in 
computer  control  systems  that  use  measurement  of  the  ball  screw  rotation  to  infer  the  positioning 
assuming  the  constant  calibrated  pitch  of  the  ball  screw.  This  type  of  metrology  system  is 
currently  used  for  strut  length  positioning  of  the  NIST  Octahedral  Hexapod  machine  tool. 

A motor  is  used  to  rotate  the  ball  screw.  The  rated  continuous  power  is,  typically,  several 
kilowatts.  A ball  screw-ball  nut  system  is  usually  90  percent  efficient.  That  is,  90  percent  of  the 
motor  power  is  transferred  from  the  motor  input  to  the  ball  nut  to  move  the  mechanical  load. 

The  ten  percent  loss  in  power  at  the  ball  screw-ball  nut  contact  is  friction  in  the  system  and  the 
fiiction  is  exhibited  as  heat  generation.  The  heat  generated  in  the  ball  nut  is  transferred  to  both 
the  machine  attachment  and  the  ball  screw.  This  heat  transfer  causes  the  ball  screw  to  expand 
and  changes  the  ball  screw  pitch.  The  change  in  pitch  results  in  errors  in  positioning.  Based 
upon  the  efficiency  of  the  system,  the  thermal  power  transferred  to  the  ball  screw  may  be  on  the 
order  of  a kilowatt. 

The  analysis  described  in  this  report  is  an  approximation  to  the  physical  situations  discussed 
above.  The  purpose  of  the  analysis  is  two  fold:  a simulation  model  to  estimate  the  temperature- 
dependent  strut  length  errors  of  the  NIST  Octahedral  Hexapod  and  ultimately  the  tool-tip 


2 


positioning  errors  of  the  machine  and  to  assist  in  the  interpretation  of  ball  screw  temperature 
measurements  on  the  operating  machine. 

THE  NIST  OCTAHEDRAL  HEXAPOD 

Figure  1 is  a solid  model  representation  of  the  NIST  Octahedral  Hexapod  machine  tool  [1].  The 
tool  spindle  is  supported  by  six  struts.  As  described  below,  each  strut  is  extended  and  retracted 
in  a coordinated  manner  to  position  the  tool  tip  for  machining  a part.  The  tool  tip  position  and 
the  orientation  of  the  tool  axis  relative  to  the  part  is  dependent  upon  the  lengths  of  all  the  struts. 
The  computer  control  of  the  machine  utilizes  data  describing  the  calibrated  geometry  of  the 
machine  tool.  This  calibrated  geometry  comprises,  in  part,  the  locations  in  space  of  each  end  of 
each  strut  and  the  length  of  each  strut  for  the  “home”  position  of  the  machine  tool.  The 
calibrated  geometry  was  established  by  the  manufacturer  upon  the  machine’s  assembly  at  NIST. 

Figure  2 illustrates  the  strut  geometry  of  the  NIST  Octahedral  Hexapod.  On  the  right  end  of  this 
figure  a large  ball  joint  is  illustrated  (329  mm  diameter’).  This  is  the  strut  end  supported  by  the 
octahedral  frame  of  the  machine  tool.  Rigidly  attached  to  the  large  ball  joint  is  the  outer  sleeve 
of  the  strut.  As  illustrated  in  Figure  3,  the  outer  sleeve  is  2062.4  mm  long.  An  inner  sleeve 
slides  within  the  outer  sleeve.  The  left  hand  end  of  the  figure  illustrates  an  80  mm  diameter  ball 
joint.  This  ball  joint  is  rigidly  attached  to  the  sliding  inner  sleeve  and  connects  the  strut  to  the 
platform  supporting  the  tool  spindle.  As  illustrated,  the  length  of  the  sliding  inner  sleeve  is 
2037.4  mm.  Also  illustrated  in  this  figure  at  the  right  hand  end  of  the  sliding  inner  sleeve  is  a 
260  mm  long  area  denoted  by  an  “x”  symbol.  This  is  the  ball  nut  and  it  is  rigidly  attached  to  the 
inner  sleeve.  The  sliding  inner  sleeve  is  keyed  to  the  outer  sleeve  so  that  the  inner  sleeve  cannot 
rotate  with  respect  to  the  outer  sleeve. 

Figure  3 illustrates  the  strut  in  three  configurations  of  strut  extension.  This  extension  is 
mechanically  realized  by  rotating  a ball  screw  which  drives  the  ball  nut.  In  Figures  2 and  3 the 
ball  screw  is  illustrated  by  the  /•/  symbol.  From  the  center  of  the  large  ball  joint  of  the  strut,  the 
ball  screw  length  is  1912.4  mm.  As  indicated  in  Figure  3 (see  notes  in  upper  left;  hand  comer), 
mechanical  limitations  result  in  a limited  stmt  stroke  of  1433  mm. 

The  stmt  ball  screw  is  driven  by  a mechanical  coupling  to  a drive  motor  at  the  right  hand  end  of 
the  stmt.  The  coupling  and  motor  are  not  illustrated  in  Figures  2 and  3.  Also,  the  location  of  the 
switches  indicating  the  stmt  “home”  position  is  not  illustrated.  These  switches  are  located  on  the 
outer  sleeve  of  the  stmt  on  the  right  hand  side  and  towards  the  large  ball  joint.  The  stmt  “home” 
position  is  a configuration  similar  to  that  illustrated  at  the  top  of  Figure  3. 

For  the  NIST  Octahedral  Hexapod,  the  ball  screw  of  each  stmt  is  manufactured  with  a 10  mm 
diameter  hole  (not  shown)  through  the  length  of  the  ball  screw.  Also,  the  clearance  hole 
continues  through  the  coupling  to  the  drive  motor  and  the  drive  motor  axis.  The  purpose  of  this 
hole  is  to  allow  for  installation  of  a laser-based  stmt  length  metrology  system  for  each  stmt. 


’all  dimensions  are  design  values 
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FIGURE  1 The  NIST  Octahedral  Hexapod  Machine  Tool 
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FIGURE  2 NIST  Octahedral  Hexapod  Strut:  Basic  Dimensions 
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Presently,  the  stmt  length  metrology  is  based  upon  a resolver  fixed  to  the  stmf  drive  motor  shaft 
and  stmt  length  is  inferred  from  the  motor  rotation  and  the  calibrated  fixed  pitch  of  the  ball 
screw. 

Friction  of  sliding  contact  is  exhibited  by  heat  generation.  Three  potentially  significant  sources 
of  heat  generation  are  evident  in  the  stmt  configuration  illustrated  in  Figures  2 and  3.  The 
sliding  contact  between  the  inner  and  outer  sleeve  of  the  stmt  results  in  heating  of  both 
components.  This  heating  causes  both  components  to  experience  a temperature  change  and 
thermal  expansion  which  is  both  distance  dependent  (along  the  stmt  component)  and  time 
dependent.  The  sliding  contact  area  between  the  inner  sleeve  and  the  outer  sleeve  is  the  first  400 
mm  measured  from  the  left  hand  end  of  the  outer  sleeve.  This  thermal  expansion  may  result  in 
two  sources  of  length  error  for  the  stmt  extension:  expansion  of  the  2037.4  mm  long  iimer  sleeve 
and  “drift”  of  the  reference  location  of  the  “home”  position  switches  fixed  to  the  outer  sleeve. 

For  the  present  configuration,  the  computer  control  of  the  NIST  Octahedral  Hexapod  is  ignorant 
of  these  effects. 

The  second  source  of  potential  sliding  contact  heat  generation  is  the  swivel  of  the  stmt  ball  joints 
within  their  mechanical  seats.  For  these  joints,  the  rotation  rates  are  rather  lower  in  magnitude 
when  compared  to  the  stmt  extension  rates.  This  aspect  of  the  heat  generation  sources  is  not 
considered  here. 

The  third  potential  source  of  heat  generation  is  the  ball  nut  friction.  The  heat  generation  within 
the  rolling  elements  of  the  ball  nut  is  transferred  to  both  the  sliding  inner  sleeve  and  the  ball 
screw.  The  heat  transferred  to  the  iimer  sleeve  will  contribute  to  its  thermal  expansion  in 
addition  to  the  sliding  contact  with  the  outer  sleeve.  The  ball  nut  heat  transferred  to  the  ball 
screw  will  result  in  thermal  expansion  of  the  ball  screw.  The  thermal  expansion  of  the  ball  screw 
results  in  a nonuniform  change  in  the  ball  screw  pitch.  The  nonuniformity  is  in  both  position 
along  the  ball  screw  and  in  time.  The  computer  control  of  the  NIST  Octahedral  Hexapod  is  also 
ignorant  of  these  effects. 

The  NIST  Octahedral  Hexapod  utilizes  excellent  machine  tool  design.  Provisions  are 
incorporated  in  the  design  to  maintain  well  lubricated  surfaces  for  all  sliding  components.  The 
lubrication  is  provided  to  minimize  the  effects  of  heat  generation  by  sliding  contact.  However 
minimal,  the  effects  of  thermal  expansion  must  be  quantified  in  order  to  characterize  the  tool-tip 
positioning  of  the  machine.  The  present  analysis  is  one  step  along  the  path  to  characterization  of 
this  machine  tool  configuration. 

HEAT  CONDUCTION/TRANSFER  MODEL 

The  heat  conduction/transfer  model  described  in  this  report  is  illustrated  in  Figure  4^.  Appendix 
A presents  the  mathematical  details  for  obtaining  analytical  expressions  for  the  temperature  , 


^Beginning  with  Figure  4,  the  figures  are  placed  at  the  end  of  the  text  so  as  not  to 
intermpt  the  commentary. 
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T(z,t),  and  the  thermal  expansion,  w(z,t),  of  the  rod.  Appendix  B presents  the  details  for 
estimating  the  forced  convection  heat  transfer  coefficient,  h.  The  discussion  here  uses  these 
results  without  explicity  trudging  through  the  mathematical  details. 

The  mechanical  system  illustrated  in  Figure  4 is  a one-dimensional  (z  is  the  dimension)  heat 
conduction  problem  [2].  The  temperature  distribution,  T(z,t),  must  satisfy  a partial  differential 
equation  subject  to  an  initial  condition  (state  of  the  system  at  time  t = 0)  and  a boundary 
condition  (the  heat  input  condition  at  z = 0).  The  initial  condition  applies  for  all  z > 0 and  the 
boundary  condition  applies  for  all  t > 0. 

Once  the  expression  for  the  temperature  is  obtained,  the  thermal  expansion,  w(z,t),  is  obtained  by 
integrating  the  expression  for  the  temperature.  This  form  of  the  thermal  expansion  is  specifically 
a “stress  free”  expansion.  That  is,  there  are  no  constraints  imposed  on  the  rod  to  oppose  the 
expansion. 

Now,  we  discuss  the  two  physical  parameters  that  characterize  both  the  temperature  and  thermal 
expansion  of  the  rod.  First,  we  consider  the  temperature  and  then  the  thermal  expansion.  To  do 
this  it  is  necessary  to  state  the  form  of  the  governing  partial  differential  equation. 


The  governing  partial  differential  equation  is  derived  in  Appendix  A and  is  given  by  Eq(A-5b). 
The  equation  is: 


d^T  . dT 

K - X(T-T  ) = — : z > 0 and  r > 0 

dz^  ” 


(1) 


For  this  equation  T(z,t)  is  the  temperature  (solution  of  the  equation)  and  is  a function  of  the  axial 
distance,  z,  and  time,  t.  The  parameters  characterizing  the  problem  are  k,  A,  and  T„  which  are 
assumed  constant.  As  indicated  in  Figure  4,  T„  is  the  constant  ambient  temperature  of  the  air 
surrounding  the  rod  and  is  the  temperature  used  to  specify  the  initial  condition  for  Eq  (1). 
Discussion  will  focus  upon  the  two  really  important  parameters:  k and  A.  These  are  dimensional 
parameters  (i.e.,  their  magnitude  depends  upon  the  physical  units  used  to  define  mass,  length, 
time  and  temperature). 

The  physical  units  used  here  are:  mass  in  kilograms  (kg),  length  in  meters  (m),  time  in  seconds 
(s),  and  temperature  in  degrees  Celsius  (°C).  Each  term  in  Eq  (1)  has  the  physical  dimension  of 
°C/s.  As  described  in  Appendix  A,  the  first  term  in  Eq(l)  (the  term  containing  k)  is  the  time  rate 
of  change  of  thermal  energy  due  to  heat  conduction  along  the  rod  at  the  location  z and  time  t of 
the  system.  The  second  term  in  Eq(l)  (the  term  containing  A)  is  the  time  rate  of  change  of 
thermal  energy  due  to  the  dissipation  of  heat  along  the  length  of  the  rod.  The  right  hand  side  of 
Eq(l)  is  the  time  rate  of  change  of  thermal  energy  at  the  location  z and  time  t of  the  system. 

The  parameter  k is  called  the  thermal  diffusivity  and  depends  upon  the  material  properties  of  the 
rod.  The  physical  units  of  k are  m^/s.  This  parameter  scales  (i.e.,  relates)  the  spatial  coordinate  z 
(m)  and  the  temporal  coordinate  t (s). 
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The  parameter  X is  called  the  cooling  coefficient  and  it  is  directly  proportionaTto  the  convection 
heat  transfer  coefficient,  h,  and  is  also  a ftmction  of  the  cross-section  shape  of  the  rod  and  the 
material  properties  of  the  rod  (see  Eq  (A-5)).  This  is  a dissipation  rate  and  has  physical 
dimensions  of  1/s.  This  parameter  scales  the  temporal  coordinate,  t.  The  parameter  X 
characterizes  the  time  constant  for  the  temperature  variation  of  the  system. 

The  parameter  / (A./k)  has  physical  dimension  1/m.  As  indicated  in  the  rather  complicated 
expressions  for  temperature  given  in  Appendix  A,  the  parameter  V (A/k)  characterizes  the  spatial 
constant  for  the  temperature  variation. 

Appendix  B describes  the  steps  necessary  to  estimate  the  parameter  X and  presents  the  ranges  of 
the  magnitudes  of  the  parameters  characterizing  the  problem.  For  a circular  cross  section  steel 
rod  cooled  by  low  velocity  air  flow  at  atmospheric  temperature  and  pressure,  the  ranges  of 
magnitude  of  k and  A.  are: 


1 T 0'^  mVs  < K < 2T  0'^  m^/s 
and 

1.2-10-'  l/s<  A <7.2-10-^  1/s. 

From  these  ranges  one  may  infer  the  range  of  the  spatial  constant  as: 

2.45  1/m  < / (A./k)  < 26.8  1/m. 

The  constants  X and  / (A,/k)  indicate  that  temperature  changes  occur  over  several  thousands  of 
seconds  (i.e.  hours)  and  over  distances  of  less  than  one  meter,  for  the  stated  conditions. 

The  thermal  expansion,  w(z,t),  is  characterized  by  the  above  parameters  and  an  additional 
parameter,  a:  the  coefficient  of  thermal  expansion.  The  physical  dimensions  of  a are  m/m°C  or 
1/°C.  For  steels  used  in  precision  machines  [3]  the  range  of  the  coefficient  of  thermal  expansion 
is: 

almost  zero  < a < IS-IO"^  1/°C. 

For  a ball  screw,  the  upper  limit  for  a is  more  typical.  Measuring  thermal  expansion  in 
micrometers  (^im)  it  is  seen  that  a few  degress  Celsius  change  in  temperature  can  result  in  a 
significant  thermal  expansion.  Since  the  thermal  expansion  involves  the  integral  of  the 
temperature,  local  changes  in  temperature  along  the  rod  result  in  larger  cumulative  displacements 
(i.e.,  position  errors)  along  the  rod. 

Example  calculations  illustrate  the  above  comments. 

EXAMPLE  CALCULATIONS 

Appendix  A presents  the  mathematical  expressions  for  the  temperature,  T(z,t),  and  thermal 
expansion,  w(z,t),  of  the  rod.  These  expressions  assume  an  initial  temperature  of  the  rod  at  the 
ambient  temperature  of  the  air  (the  initial  condition  for  Eq(l)).  Two  physical  forms  of  the  heat 
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input  to  the  rod  are  considered:  constant  temperature  input  and  constant  heat  flux  input.  These 
are  the  two  boundary  conditions  for  Eq(l).  Hence,  two  problems  are  considered:  constant 
temperature  input  and  constant  heat  flux  input.  This  form  of  input  is  denoted  mathematically  as 
a step  function  input.  The  resulting  expressions  for  temperature  and  thermal  expansion  of  the  rod 
are  the  step  response  functions  of  the  system.  These  step  response  functions  may  be  used  to 
obtain  temperature  and  thermal  expansion  for  more  complicated  heat  input  conditions  than  those 
considered  in  this  report  [4]. 

The  actual  heat  input  to  a ball  screw  in  contact  with  a moving  ball  nut  is  (as  an  understatement)  a 
very  complicated  problem.  In  the  present  discussion,  the  two  physical  forms  of  heat  input  are 
considered  since  the  actual  input  may  never  be  either  measurable  or  known.  What  can  be 
measured  are  temperatures  at  a distance  away  from  the  input  and  overall  thermal  expansion 
(position  errors)  of  the  system.  Hence,  the  results  in  this  report  can,  at  best,  assist  in  the 
interpretation  of  experimental  data. 

PARAMETERS  FOR  EXAMPLE  CALCULATIONS 

As  described  above,  the  two  parameters  defining  the  spatial  and  temporal  characteristics  are  k 
and  X.  To  establish  the  amplitude  of  the  system  step  response  it  is  necessary  to  specify  the  initial 
ambient  temperature  and  the  magnitude  of  the  heat  input. 

For  illustration,  the  rod  is  assumed  to  be  303  stainless  steel  at  an  initial  temperature  of  23  °C  and 
atmospheric  pressure.  The  diameter  of  the  rod  is  80  mm  and  the  cooing  air  velocity  is  0.1  m/s. 
For  the  thermal  expansion,  it  is  assumed  that  the  rod  scale  (i.e.,  ball  screw  pitch)  is  calibrated  at 
20 °C.  Hence,  the  rod  in  its  initial  state  of  23  °C  exhibits  a thermal  expansion  due  to  the  3 °C 
temperature  change  between  the  scale  calibration  temperature  and  the  ambient  temperature. 

The  data  for  the  calculations  are  [3]: 


K = 4.1T0'^  mVs 
A = 5.0-10-'  1/s 
a = 17.2-10-^  \rc 
ambient  temperature  = 23  °C 
calibration  temperature  = 20°C 
step  temperature  input  = 20°C 
step  heat  flux  input  = 1 13 1 .6  W/m^ 

The  value  for  the  step  heat  flux  input  is  specified  so  that  the  asymptotic  temperature  rise  at  z = 0 
for  constant  heat  flux  input  is  20°C.  (See  Eq(A-14)  and  the  limiting  expressions  for  the 
temperature  rise.)  Also,  the  thermal  conductivity  of  303  stainless  steel,  k,  is  16.2  W/m°C.  The 
magnitude  of  the  step  response  inputs  are  taken  only  for  illustration  and  result  in  a final  (or 
asymptotic)  temperature  of  43  °C. 
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STEP  RESPONSE  FOR  CONSTANT  TEMPERATURE  INPUT 


The  temperature  is  a function  of  distance  along  the  rod  and  time.  Figure  5 is  a three  dimensional 
plot  of  the  temperature  as  a function  of  distance  along  the  rod  in  meters  and  time  in  minutes.  The 
distance  scale  covers  the  range  0 m to  0.250  m.  The  time  scale  covers  the  range  0 minutes  to  480 
minutes  (8  hours).  The  constant  temperature  of  23  °C  along  the  plane  t = 0 minutes  is  the  initial 
temperature  of  the  rod.  Considering  the  plane  z = 0,  it  is  seen  that  the  constant  temperature  input 
of  20°C  is  realized  immediately  and  remains  constant  as  time  increases.  The  final  temperature  of 
the  system  at  z = 0 is  (23  °C  + 20°C)  43 °C. 

Considering  the  plane  t = 480  minutes  (8  hours),  the  temperature  distribution  along  the  length  of 
the  rod  exhibits  an  exponential  decay.  (See  the  limiting  expressions  for  temperature  following 
Eq(A-l  1).)  The  spatial  decay  rate  for  this  example  is  / (A/k)  = 3.49  1/m.  That  is,  after  steady 
temperatures  are  realized  (very  large  times),  at  z = 1/3.49  = 0.286  m the  temperature  has 
decreased  by  a factor  of  e'’  = 0.37  or  AT  = 0.37(20)  = 7.4 °C  above  ambient  temperature. 

Considering  the  plane  z = 0.250m,  the  temperature  rise  exhibits  a smooth  exponential-like 
increase  once  significant  thermal  energy  reaches  this  location.  The  explicit  functional  form  is 
given  by  Eq(A-l  1).  The  time  constant  for  this  example  is  \/X  = 20,000  seconds  (333.3  minutes). 

Figure  6 presents  contours  of  constant  temperature  as  a function  of  distance  along  the  rod  and 
time.  This  figure  indicates  that  the  rod  never  realizes  a constant  temperature  distribution  along 
its  length.  This  is  expected  since  the  rod  is  cooled  along  its  length. 

Figure  7 is  the  step  response  of  the  thermal  expansion  for  a constant  temperature  input.  The 
thermal  expansion  is  measured  relative  to  the  plane  z = 0 and  is,  by  definition,  zero  at  z = 0.  The 
plane  t = 0 is  the  thermal  expansion  at  an  instant  before  the  step  temperature  input  is  applied. 

The  constant  increase  in  thermal  expansion  is  due  to  the  3°C  temperature  difference  between  the 
ambient  temperature  of  the  system  and  the  assumed  calibration  temperature  of  the  rod.  (For 
computational  reasons,  the  initial  time  illustrated  in  Figure  7 is  0.1  minutes.) 

Figure  8 is  a contour  plot  of  constant  thermal  expansion  as  a function  of  distance  along  the  rod 
and  time.  This  figure  illustrates  that  after  480  minutes  of  heating  the  first  200  mm  of  rod  are 
approaching  steady  (time  independent)  conditions.  Beyond  200  mm  of  rod  length,  the  thermal 
expansion  is  slowly  changing.  For  the  first  60  minutes  of  heating.  Figure  8 indicates  that  the 
thermal  expansion  is  dominated  by  the  3 °C  temperature  difference  between  the  ambient 
temperature  and  the  calibration  temperature.  This  illustrates  the  importance  of  calibrating  a 
machine  tool  in  the  ambient  environment  that  it  will  be  operated  or  making  real-time  correction 
for  this  error. 

STEP  RESPONSE  FOR  CONSTANT  HEAT  FLUX  INPUT 

As  described  in  Appendix  A,  the  heat  flux  is  a function  of  distance  along  the  rod  and  time.  The 
heat  flux  satisfies  the  same  partial  differential  equation  as  does  the  temperature,  Eq(l).  For  the 
rod  initially  at  the  uniform  ambient  temperature,  the  initial  heat  flux  is  zero.  The  step  increase  of 
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a constant  heat  flux  input  is  mathematically  identical  to  the  constant  temperature  input  discussed 
above.  Hence,  the  functional  form  for  the  heat  flux  step  response  is  identical  to  the  temperature 
step  response  relative  to  the  constant  ambient  temperature. 

Figure  9 is  a three  dimensional  plot  of  heat  flux  response  for  the  constant  heat  flux  input  of 
1131.6  W/m^.  The  magnitude  of  the  input  is  selected  for  illustration  so  that  the  reader  may 
easily  compare  the  step  temperature  response  and  the  step  thermal  expansion  response  for  heat 
flux  input  to  the  results  (Figures  5 through  8)  for  constant  temperature  input.  Comparing  Figures 
5 and  9,  it  is  seen  that  the  shape  of  the  surfaces  is  the  same.  All  of  the  discussion  for  Figure  5 
concerning  time  constants,  spatial  decay  and  limiting  conditions  are  applicable  to  Figure  9. 

Figure  10  is  a constant  contour  plot  of  heat  flux  as  a function  of  distance  along  the  rod  and  time. 
This  result  is  comparable  to  Figure  6 for  the  temperature  step  response. 

For  the  step  heat  flux  input,  the  temperature  response  is  obtained  by  integrating  the  heat  flux  with 
respect  to  the  spatial  coordinate  z (see  Eq  (A- 13)).  This  step  is  necessary  since  the  heat  flux  is 
proportional  to  the  spatial  temperature  gradient.  The  proprotionality  constant  is  the  thermal 
conductivity,  k,  of  the  rod  material.  For  this  example,  k = 16.2  W/m°C. 

Figure  11  is  the  temperature  response  of  the  rod  due  to  the  step  heat  flux  input.  This  result 
should  be  compared  to  Figure  5 for  the  step  temperature  input.  For  the  plane  z = 0,  the  step  heat 
flux  input  causes  a gradual  temperature  rise  from  the  ambient  temperature.  For  the  plane  z = 
0.250  m,  the  temperature  rise  is  also  gradual  and  is  similar  to  that  illustrated  in  Figure  5.  The 
explicit  functional  form  is  given  by  Eq(A-14).  For  the  plane  t = 480  minutes,  the  temperature 
response  is  a spatial  decay  along  the  length  of  the  rod  and  is  functionally  identical  to  that 
discussed  for  Figure  5.  That  is,  after  the  steady  temperature  condition  has  been  realized  at  z = 0, 
the  spatial  temperature  distribution  along  the  rod  is  identical  for  both  step  temperature  input  and 
step  heat  flux  input. 

Figure  12  presents  contours  of  constant  temperature  response  due  to  the  step  heat  flux  input. 

This  result  should  be  compared  with  Figure  6 for  the  step  temperature  input.  Hence,  it  is  seen 
that  the  temperature  change  along  the  rod  for  any  time  is  more  gradual  than  that  for  the  step 
temperature  input.  Although  the  data  for  this  example  was  selected  so  that  the  temperature  rise 
would  asymptotically  reach  43  °C  at  z = 0,  it  is  seen  in  Figure  12  that  after  480  minutes  the 
temperature  response  has  not  reached  its  limiting  value. 

Figure  13  is  the  thermal  expansion  response  of  the  rod  due  to  the  step  heat  flux  input.  This  result 
should  be  compared  to  Figure  7 for  the  step  temperature  input.  The  comparison  illustrates  that 
the  thermal  expansion  is  more  gradual  for  the  step  heat  flux  input  as  compared  to  the  step 
temperature  input.  This  is  to  be  expected  since  the  temperature  response  is  more  gradual  when 
compared  to  the  step  temperature  input.  The  increasing  thermal  expansion  at  time  t = 0 results 
from  the  3 °C  temperature  difference  between  the  23  °C  ambient  temperature  and  the  20 °C  rod 
calibration  temperature  assumed  for  this  example. 
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Figure  14  is  a plot  of  contours  of  constant  thermal  expansion  as  a function  of  distance  along  the 
rod  and  time  due  to  the  step  heat  flux  input.  This  result  should  be  compared  with  Figure  8 for  the 
step  temperature  input.  For  the  step  heat  flux  input,  the  contours  of  constant  thermal  expansion 
are  more  gradual  both  in  distance  along  the  rod  and  in  time  as  compared  to  the  step  temperature 
input.  As  indicated  in  Figure  14,  after  480  minutes  of  heating,  steady  state  conditions  are  still 
not  realized. 

LIMITING  COOLING  CONDITIONS 

The  cooling  parameter,  X,  is  the  time  rate  of  heat  dissipation  per  unit  length  of  the  rod.  For  a 
circular  cross  section  steel  rod  cooled  by  low  velocity  air  flow  at  atmospheric  temperature  and 
pressure,  the  numerical  range  for  X is: 

1.2*10-^  1/s  <X<  7.2*10-^  1/s 

It  is  instructive  to  calculate  the  thermal  expansion  corresponding  to  the  extreme  values  of  this 
range.  For  this  example  the  data  are  taken  as  described  previously  for  constant  heat  flux  input 
(Figures  9 through  14).  The  only  change  is  the  value  of  X.  Two  values  of  X are  discussed: 

Low  Cooling:  X,o^=  1.2*10'^  1/s 

and 

High  Cooling:  X^igh  = 7.4T0'^  1/s. 

For  each  of  these  cases,  two  results  are  discussed:  the  temperature  response  of  the  rod  and  the 
thermal  expansion  of  the  rod.  These  results  should  given  the  reader  an  appreciation  of  the 
possible  significance  of  cooling  effects  as  estimated  by  the  model.  For  the  results  presented  in 
Figures  9-14,  the  cooling  rate  is  X = 5.0*  10'^  1/s. 

For  low  cooling,  Figure  15  is  the  plot  of  contours  of  constant  temperature  for  distance  along  the 
rod  and  time.  Figure  15  should  be  compared  to  Figure  12.  Two  differences  are  evident  between 
Figure  12  and  Figure  15:  The  temperature  distribution  along  the  rod  at  fixed  time  and  the 
asymptotic  temperature  at  the  rod  end  (z  = 0).  Figure  15  corresponds  to  a lower  cooling  rate  as 
compared  to  Figure  12.  Hence,  at  any  time  following  initiation  of  heat  flux  input,  the 
temperature  distribution  along  the  rod  length  is  larger  in  magnitude  in  Figure  15  than  in  Figure 
12.  For  Figure  15,  the  asymptotic  temperature  exceeds  45  °C  (for  Figure  12,  the  heat  flux  input 
was  selected  to  result  in  a 20°  C temperature  rise  above  ambient).  The  asymptotic  temperature 
rise  above  ambient  is  given  by  the  expression  for  T*  in  Eq(A-14e).  For  the  cooling  rate  1.2*10'^ 
of  Figure  15,  the  asymptotic  temperature  rise  is  40.8  °C  above  ambient  (63.8  °C  total  for  an 
ambient  temperature  of  23  °C).  This  illustrates  that  for  a constant  heat  flux  input  and  low 
cooling,  the  temperature  distribution  along  the  rod  is  distinctly  different  from  a more  “typical” 
condition  as  illustrated  in  Figure  12. 
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Figure  16  is  the  plot  of  contours  of  constant  thermal  expansion  corresponding  to  the  “low” 
cooling  temperature  distribution  of  Figure  15.  The  thermal  expansion  contours  of  Figure  16 
should  be  compared  with  the  contours  of  Figure  14.  In  this  comparison,  the  reader  should 
recognize  that  the  thermal  expansion  is  more  severe  for  the  for  the  lower  cooling  rate  illustrated 
in  Figure  16. 

Now,  we  consider  a high  cooling  rate  corresponding  to  A = 7.4T0‘^  1/s.  This  rate  corresponds  to 
a much  higher  velocity  forced  air  cooling  than  previously  described.  The  feasibility  of  achieving 
such  a high  cooling  rate  is,  perhaps,  unrealistic.  However,  the  purpose  here  is  to  illustrate 
extreme  effects. 

Figure  17  is  the  plot  of  contours  of  constant  temperature  for  A = 7.4T0’^  1/s.  This  plot  should  be 
compared  to  Figures  12  and  15.  Figure  17  illustrates  a very  dramatic  difference  in  temperature 
response.  First,  for  this  “high”  cooling  example,  steady  state  temperature  distribution  is  achieved 
within  120  minutes  of  initiation  of  heat  flux  input,  and  the  ambient  (23  °C)  temperature  is 
realized  for  rod  locations  within  approximately  0.3  m from  the  heat  input.  Additionally,  the 
asymptotic  temperature  rise  for  the  rod  is  5.3  °C  above  ambient  rather  than  the  20  °C  rise  for 
Figure  12  and  the  40.8°C  for  Figure  15.  For  each  of  these  results,  the  reader  is  reminded  that  the 
input  heat  flux  is  constant  at  1 13 1 .6  W/m^. 

Figure  18  is  the  plot  of  contours  of  constant  thermal  expansion  for  the  temperature  response  of 
Figure  17.  The  response  illustrated  in  Figure  18  should  be  compared  with  the  results  of  Figures 
14  and  16.  The  “high”  cooling  response  illustrated  in  Figure  18  is  also  dramatically  different 
from  the  thermal  expansion  response  illustrated  in  Figures  14  and  16.  Figure  18  illustrates  that 
steady  state  conditions  are  realized  within  the  first  1 20  minutes  of  initiating  the  constant  heat  flux 
input.  Further,  for  a fixed  location  along  the  rod,  the  thermal  expansion  for  “high”  cooling 
(Figure  18)  is  approximately  1/4  the  thermal  expansion  for  “low”  cooling  (Figure  16). 

In  summary.  Figures  15  and  16  illustrate  temperature  and  thermal  expansion  response, 
respectively,  for  “low”  cooling  along  the  rod  length.  Figures  17  and  18  are  the  corresponding 
response  characteristics  for  “high”  cooling.  For  both  cases,  the  illustrated  response  is  a result  of 
a step  heat  flux  input  to  the  rod  of  constant  magnitude. 

REALITY  CHECK 

The  above  results  are  illustrations  of  an  idealized  system  response  to  two  different  input 
functions.  These  results  are  based  upon  calculations  utilizing  complicated  mathematical 
expressions.  The  question  as  to  the  reality  of  either  input  function  remains  open.  Physically,  it  is 
the  author’s  opinion  that  the  step  heat  flux  model  is  closer  to  reality  than  the  step  temperature 
input  model. 

This  opinion  is  based  upon  energy  considerations.  To  physically  realize  a constant  temperature 
input,  one  must  initially  supply  an  infinite  amount  of  power  or  energy  to  the  system  and 
gradually  decrease  this  input  to  match  the  energy  dissipated  to  cooling  while  maintaining  the 
temperature  at  z = 0 constant.  This  is  impossible. 
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A constant  heat  flux  input  implies  a constant  power  or  energy  source.  How  this  may  be  realized 
mechanically  or  measured  in  reality  is  also  an  open  question.  However,  since  the  implication  is  a 
finite  heat  source,  the  author  concludes  that  the  step  heat  flux  input  is  a more  realistic  condition. 

Mathematically,  knowing  the  step  input  response  of  a system  allows  one  to  compute  the  response 
to  arbitrary  time  dependent  inputs  [4].  In  this  sense,  either  model  may  be  appropriate  although 
the  author  still  learns  towards  the  step  heat  flux  input  model.  Response  to  an  arbitrary  time 
dependent  input  apparently  requires  the  numerical  integration  of  the  complicated  mathematical 
expressions  given  in  Appendix  A.  Either  approach  requires  the  reality  of  metrology  of  the 
system  response. 

METROLOGY  CONSIDERATIONS 

As  mentioned  above,  an  infinite  supply  of  energy  is  impossible.  This  implies  that  a constant  step 
temperature  input  condition  is  impossible.  From  a metrology  standpoint,  the  measurement  of 
heat  flux  for  mechanical  contact  of  two  surfaces  is  only  slightly  less  than  impossible. 
Measurement  of  temperature,  especially  under  time-varying  conditions,  is  rather  easily 
accomplished. 

The  results  described  in  this  report  are  a guide  to  the  conduct  of  metrology  experiments  and  the 
interpretation  of  the  experimental  data.  On  the  negative  side  one  may  not  be  aware  of  details  of 
the  heat  input.  On  the  positive  side  one  has  at  least  a few  minutes  to  acquire  metrology  data  with 
the  knowledge  that  the  conditions  are  not  significantly  changing  with  time. 
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FIGURE  4 Heat  Conduction/Transfer  Model 


CONSTANT  TEMPERATURE  INPUT 


FIGURE  5 Step  Temperature  Response:  20°C  Step  Input  and  23 °C 

Ambient  Temperature 
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time,  minutes 


CONTOURS  OF  CONSTANT  TEMPERATURE,  deg  C 


FIGURE  6 Contours  of  Constant  Temperature  for  Step 

Temperature  Response:  20°C  Step  Input  and  23 °C 
Ambient  Temperature 
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thermal  expansion,  urn 


CONSTANT  TEMPERATURE  INPUT 


distance  along  rod,  m 


FIGURE  7 Thermal  Expansion  for  a Step  Temperature  Input: 

20  °C  Step  Input,  23  °C  Ambient  Temperature  and  20°  C 
Rod  Calibration  Temperature 
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CONTOURS  OF  CONSTANT  THERMAL  EXPANSION,  um 


FIGURE  8 Contours  of  Constant  Thermal  Expansion  for  a Step 

Temperature  Input:  20  °C  Step  Input,  23  °C  Ambient 
Temperature  and  20  °C  Rod  Calibration  Temperature 
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heat  flux,  W/m'^2 


CONSTANT  HEAT  FLUX  INPUT 


distance  along  rod,  m 
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FIGURE  9 Step  Heat  Flux  Response  (Heat  Flux  Step  Amplitude  to 

Yield  20 °C  Asymptotic  Temperature  Rise) 
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CONTOURS  OF  CONSTANT  HEAT  FLUX,  W/ni''2 


FIGURE  10  Contours  of  Constant  Heat  Flux  for  Step  Heat  Flux 

Input 
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CONSTANT  HEAT  FLUX  INPUT 


FIGURE  1 1 Temperature  Response  to  Step  Heat  Flux  Input: 

Asymptotic  20  °C  Temperature  Rise  and  23  °C  Ambient 
Temperature 
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time,  minutes 


CONTOURS  OF  CONSTANT  TEMPERATURE,  deg  C 


distance  along  rod,  m 


FIGURE  1 2 Contours  of  Constant  Temperature  for  a Step  Heat  Flux 

Input:  Asymptotic  20°C  Temperature  Rise  and  23 °C 
Ambient  Temperature 
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thermal  expansion,  urn 


CONSTANT  HEAT  FLUX  INPUT 


distance  along  rod,  m 


FIGURE  13  Thermal  Expansion  for  a Step  Heat  Flux  Input; 

Asymptotic  20  °C  Temperature  Rise,  23  °C  Ambient 
Temperature,  and  20  °C  Rod  Calibration  Temperature 
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time,  minutes 


CONTOURS  OF  CONSTANT  THERMAL  EXPANSION,  um 


FIGURE  14  Contours  of  Constant  Thermal  Expansion  for  a Step 

Heat  Flux  Input:  Asymptotic  20°C  Temperature  Rise, 
23  °C  Ambient  Temperature,  and  20  °C  Rod  Calibration 
Temperature 
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time,  minutes 


CONTOURS  OF  CONSTANT  TEM PERATURE,  deg  C 


distance  along  rod,  m 


FIGURE  1 5 Contours  of  Constant  Temperature  Response  for  “Low” 

Cooling  (A,  = 1.2T0'^  1/s)  Due  to  a Constant  Step  Heat 
Flux  Input  (Compare  to  Figures  12  and  17) 
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time,  minutes 


CONTOURS  OF  CONSTANT  THERMAL  EXPANSION,  urn 


FIGURE  1 6 Contours  of  Constant  Thermal  Expansion  Response  for 

“Low”  Cooling  (X  = 1.2T0'’  1/s)  Due  to  a Constant 
Step  Heat  Flux  Input  (Compare  to  Figures  14  and  18) 
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time,  minutes 


CONTOURS  OF  CONSTANT  TEMPERATURE,  deg  C 


distance  along  rod,  m 


FIGURE  1 7 Contours  of  Constant  Temperature  Response  for 

“High”  Cooling  (A  = 7.4*10’^  1/s)  Due  to  a Constant 
Step  Heat  Flux  Input  (Compare  to  Figures  12  and  15) 
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time,  minutes 


CONTOURS  OF  CONSTANT  THERMAL  EXPANSION,  urn 


FIGURE  1 8 Contours  of  Constant  Thermal  Expansion  Response  for 

“High”  Cooling  (A  = 7.4T0'^  1/s)  Due  to  a Constant 
Step  Heat  Flux  Input  (Compare  to  Figures  14  and  16) 


30 


REFERENCES 


1.  Falco,  J.A.  (personal  communication)  The  solid  model  illustrated  is  one  result  of  the 
simulation  modeling  of  the  NIST  Octahedral  Hexapod  Machine  Tool  by  The  Intelligent 
Systems  Division,  Manufacturing  Engineering  Laboratory,  NIST 

2.  Carslaw,  H.S.,  and  Jaeger,  J.C.:  Conduction  of  Heat  in  Solids  (2nd  Ed),  Clarendon  Press, 
Oxford  (1986)  (Chapter  4,  sections  4. 1-4.3). 

3.  Slocum,  A.H.:  Precision  Machine  Design,  Prentice  Hall  (1992)  (Table  7.3.1,  p.  335). 

4.  Boley,  B.A.  and  Weiner,  J.H.:  Theory  of  Thermal  Stresses,  John  Wiley  & Sons,  Inc. 
(1960)  (section  7.8(b)  pp.  232-234). 


31 


'>m‘m  ji 


APPENDIX  A 


MATHEMATICAL  DETAILS  OF  HEAT  CONDUCTFON  AND 
THERMAL  EXPANSION  OF  A SEMI-INFINITE  ROD 


This  appendix  presents  an  outline  of  the  derivation  of  the  differential  equation  for  heat 
conduction  in  a semi-infinite  length  rod,  the  solution  of  this  equation  for  the  initial  conditions 
and  boundary  conditions  described  below,  and  the  expressions  for  the  stress-free  thermal 
expansion  of  the  rod.  The  heat  conduction  problems  presented  here  are  completely  described  in 
Chapter  4 of  Carslaw  and  Jaeger  [A-1].  The  expressions  for  the  stress-free  thermal  expansion  of 
the  rod  are  obtained  by  integrating  the  stress-displacement  relationship  for  a rod  [A-2].  This 
approach  assumes  that  the  temperature  variation  of  the  rod  is  independent  of  the  strain  condition 
of  the  rod.  That  is,  the  temperature  variation  is  uncoupled  from  the  thermal  strain  condition. 

The  notation  used  in  this  report  for  certain  parameters  is  slightly  different  from  the  notation  used 
by  Carslaw  and  Jaeger  [A-1].  This,  however,  should  not  cause  any  confusion. 

A.  1 Differential  Equation  for  Heat  Conduction 

The  problem  considered  here  is  a long  thin  rod  exposed  to  an  ambient  temperature,  T„,  which  is 
assumed  constant.  The  rod  has  a constant  cross-section  shape  and  a varying  temperature 
distribution  T(z,t)  where  z is  the  distance  along  the  rod  and  t is  time. 

Figure  A-1  defines  the  problem  being  considered.  The  differential  equation  is  obtained  by 
considering  the  power  balance  of  the  heat  flux  across  the  surfaces  of  an  elementary  volume  of  the 
rod  material  and  the  rate  of  gain  of  heat  within  the  elementary  volume. 

The  heat  flux,  q(z,t)  per  unit  cross  sectional  area,  along  the  length  of  the  rod  is  proportional  to  the 
local  spatial  gradient  of  the  temperature  field.  The  constant  of  proportionality  is  the  thermal 
conductivity  of  the  material,  k.  The  mathematical  expression  for  the  the  flux  is: 

dT 

q{z,t)  = -k-—  {A-\) 

oz 


which  has  physical  units  of  power  per  unit  area. 

As  indicated  in  Figure  A-1,  the  net  power  gain  due  to  the  heat  flux  across  the  rod  cross-section 
areas  is: 

+k Sdz,  power  gain  due  to  heat  conduction  along  length  of  rod  (A  -2) 


A-1 


input  ^ Ambient  Temperature:  T^,  constant 


N 


CD 


A-2 


FIGURE  A- 1 Definition  of  the  Problem 


Due  to  the  temperature  difference  between  the  surface  of  the  rod  and  the  ambient  temperature, 
heat  is  lost  to  the  ambient  environment  at  the  rate 

-h(T-TJpdz,  power  lost  to  environment  (A-3) 

where  the  surface  area  is  the  perimeter  of  the  rod  cross-section,  p,  times  the  length  of  the 
element,  dz.  The  parameter  h is  discussed  in  Appendix  B. 

The  rate  of  gain  of  power  (heat  flux)  within  the  element  is  proportional  to  the  time  rate  of  change 
of  the  temperature  of  the  elementary  volume  of  the  rod.  The  rate  of  gain  of  power  is 

dT 

pc Sdz^  rate  of  gain  of  power  within  element  (-^“4) 

dt 

where  p is  the  material  density  of  the  rod  and  c is  the  specific  heat  (at  constant  pressure)  of  the 
rod  material. 

The  differential  equation  is  an  expression  of  the  equality  of  power  gain  plus  power  lost  to  the  rate 
of  gain  of  power  of  the  elementary  volume.  From  Eqs.  (A-2,  -3,  -4)  one  obtains  the  expression 
for  power  balance  as: 

(^/pc)  -(hp/pcS)(T-rj=^  (A  -5a) 

dz^  dt 


or  the  differential  equation  is  expressed  as: 

K— -X(r-rj=-^;  z > 0,  r > 0 iA-Sb) 

dz^  dt 


where  k = k/pc  thermal  diffusivity  of  the  rod  material 
X = hp/pcS  cooling  coefficient 

One  final  simplification  of  the  differential  equation  is  convenient.  We  change  the  temperature 
scale  so  that 

T(z,t)  = T*(z,t)  + T„  (A-6) 


Since  T„  is  constant,  the  differential  equation  is 

P^'p  * (37’  * 

K -XT*= ; z > 0,  r > 0 {A-1) 

dz^  dt 
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where  T*(z,t)  is  the  temperature  of  the  rod  relative  to  the  constant  ambient  temperature,  T„. 

For  no  surface  cooling,  A=0  and  equation  (A- 7)  reduces  to  the  one  dimension  heat  conduction 
equation. 


A.2  Initial  Condition  and  Boundary  Condition 

For  a solution  to  equation  (A-7)  one  must  specify  an  initial  condition  and  a boundary  condition 
on  the  function  T*. 

The  initial  condition  specifies  the  value  of  T*  at  time,  t = 0 (really  it  is  limit  T*(z,t)). 

t-  0" 

The  boundary  condition  specifies  the  value  of  T*  at  the  surfaces  of  the  rod.  For  the  present 
problem,  the  only  surface  of  the  semi-inifmite  rod  for  which  one  can  specify  the  temperature  T* 
is  the  rod  cross-section  at  z = 0 (really  it  is  limit  T*(z,t)). 

z-  0" 


Carslaw  and  Jaeger  [A-1]  is  a very  complete  compilation  of  solutions  of  heat  conduction 
problems.  It  is  a very  safe  bet  to  say  that  if  a problem  has  an  analytical  solution,  the  problem  is 
presented  in  Carslaw  and  Jaeger.  The  solutions  discussed  below  are  presented  in  Carslaw  and 
Jaeger  (Chapter  IV,  pp.  133-136). 

For  the  problem  discussed  in  the  report,  the  initial  temperature  condition  of  the  rod  is: 

limit  T(z,t)  = T„,  a constant;  z > 0 (A-8a) 

t-0" 


or 

limit  T*(z,t)  = 0;  z > 0 (A-8b) 

z-0" 


That  is,  the  rod  is  placed  in  the  environment  at  T„  and  reaches  temperature  equilibrium. 

The  boundary  condition  for  the  problem  discussed  in  the  report  is  a constant  temperature  at  the 
rod  cross-section  z = 0.  That  is,  a heat  source  is  brought  into  contact  with  the  rod  at  z = 0 and 
maintains  the  temperature  for  all  time  t > 0.  The  mathematical  expression  for  this  condition  is 

limit  T(z,t)  = To,  a constant;  z > 0 (A-9a) 

z-0" 


A-4 


or 


limit  T*(z,t)  = To  - T„=T*o,  a constant,  t > 0 (A-9b) 

z-0" 


A.3  Semi-Infinite  Rod,  Initial  Temperature  Zero,  Surface  z = 0 Maintained  at  Constant 
Temperature  for  t > 0 


The  mathematical  statement  of  this  problem  is: 


K 

dz^ 


-XT* 


dT* 

dt 


{A-\0a) 


initial  condition:  limit  T*(z,t)  = 0;  z > 0 

t-0" 

boundary  condition:  limit  T*(z,t)=T*o;  t > 0 

z-0" 


(A- 10b) 

(A- 10c) 


The  solution  of  this  problem  is  given  in  Carslaw  and  Jaeger,  Chapter  IV,  Eq.  (4.2(7))  page  135. 
The  result  is 


T*(z,t)  = T*o^(z,t) 

O(z,0  = cosh{z^|h^)-\l2e’^'^'erfizl2y|K^ 

-\l2e  ~^'^'erj{zf2yjKt-^JXt) 


(A- 11  a) 
(Allb) 


where 


dT\  is  the  error  function 


er/c(f)  = l -ery(C)  is  the  complementary  error  function. 


The  error  function  is  obviously  related  to  the  exponential  function.  It  is  useful  to  state  some 
specific  values  and  relationships  for  the  error  function: 


erf(0)  = 0 
erf(«)  = 1 
erf(-()  = -erf(0. 
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If  the  reader  is  interested,  a comprehensive  description  of  the  error  function  is  given  in  Chapter  7 
of  the  NBS  Handbook  of  Mathematical  Functions  [A-3]. 

Limiting  values  for  the  function  ^(z,t)  given  in  Eq(A-l  lb)  are: 


limit  ^ (z,t)  = 1 
z-0" 

or 

limit  T*(z,t)  = T*o 
z-0" 

(A- 11c) 

limit  ^(z,t)  = 0 

Z-oo 

or 

limit  T*(z,t)  = 0 

Z-oo 

(A- lid) 

limit  <l>(z,t)  = 0 
t-0" 

or 

limit  T*(z,t)  = 0 
t-0" 

(A- lie) 

limit  0(z,t)  = 

t-oo 

or 

limit  T*(z,t)  = T*o 

t-oo 

(A-llf) 

The  last  limit  condition  is  the  steady  or  time-independent  temperature  distribution  of  the  rod. 

A.4  Semi-Infinite  Rod,  Initial  Heat  Flux  Zero,  Surface  z = 0 Maintained  at  Constant  Heat 
Flux  for  t > 0 

Mathematically,  this  problem  is  identical  to  the  problem  described  in  A.3  above  --  only  the 
nomenclature  is  changed.  The  relation  between  the  heat  flux  and  temperature  is  given  by  Eq(A- 
1).  Clearly,  the  heat  flux,  q(z,t),  must  satisfy  the  same  differential  equation  as  the  temperature. 
The  initial  condition  and  the  boundary  condition  for  q(z,t)  are  identical  to  the  temperature 
conditions  described  in  Section  A.3.  Hence,  the  expression  for  the  heat  flux  in  the  semi-infinite 
rod  is 

q(z,t)  = q<,  $(z,t)  (A- 12) 

where  <l>(z,t)  is  given  by  Eq(A-l  lb). 

Now,  one  obtains  the  expression  for  the  temperature  T*(z,t)  by  a direct  integration  of  Eq(A-l). 
The  formal  expression  is 

fipL  dr]  = r*(ri,t)l~=-l/if~g(r),t)‘ir]  (J-I3a) 

^ z GTj  J z 

or 


{A-Ub) 
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since  limit  T*(z,t)  = 0 on  the  physical  basis  that  the  rod  is  so  long  the  locations  very  distant  from 
z 

z = 0 never  experience  a temperature  change  (i.e.,  all  the  heat  is  dissipated  due  to  the  cooling  of 
the  rod). 

The  integral  appearing  in  Eq(A-13b)  can  be  directly  obtained  in  anal>1:ical  form.  The  result  is: 
T*{z,t)  = -\jYjX{q  Jk){sir[h.(z\JXl'K) 

+ l/2e  (A-\4a) 


Limiting  values  for  the  temperature  T*(z,t)  are 

limit  T*(z,t)  = V (k/A)  (qTk)  erf(V’ (^t)) 
z^O" 

(A- 14b) 

limit  T*(z,t)  = 0 

Z-oo 

(A- 14c) 

limit  T*(z,t)  = 0 
t-0" 

(A-14d) 

limit  T*(z,t)  = / (k/X)  (qo/k)e'^’^^^’^^ 

t-oo 

(A14-e) 

The  first  limit  above  is  the  time  variation  of  the  temperature  at  z = 0 due  to  the  application  of  the 
constant  heat  flux,  q^,  at  z = 0.  For  large  times  following  initiation  of  the  heating,  the 
temperature  reaches  the  asymptotic  limit  / [K/X](qo/k).  The  last  limit  above  is  the  steady  state 
(i.e.  time  independent)  temperature  distribution  of  the  rod. 

The  reader  should  compare  the  temperature  expressions  given  in  Eqs(A-l  1)  and  Eqs(A-14).  The 
detail  expressions  are  different  — which  is  to  be  expected.  However,  the  limiting  values  of  the 
temperature  expressions  for  a large  time  are  functionally  identical.  Physically,  this  means  that 
for  large  time  after  initiating  the  heating,  the  temperature  distribution  along  the  rod  has  the  same 
functional  form  for  both  heating  conditions  considered  here.  Knowledge  of  the  steady-state 
temperature  distribution  does  not  indicate  the  nature  of  the  heat  source. 

A. 5 Stress-Free  Thermal  Expansion 
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A. 5.1  Basic  Considerations 


The  temperature  variation  along  the  length  of  the  rod  causes  a thermal  expansion  of  the  rod. 

Since  the  rod  is  long  and  slender,  the  thermal  expansion  is  one  dimensional  in  the  direction  of  the 
rod  axis.  Further,  it  is  assumed  that  the  rod  expansion  is  unconstrained  so  that  no  elastic  stress  is 
induced  in  the  rod. 

Using  the  stress-strain  relations  for  the  rod  material  and  the  strain-displacement  relation 
considering  a temperature  change  of  T(z,t),  the  thermal  expansion,  w(z,t),  is  related  to  the 
temperature  change  as: 

^ = ar(z,0  (A -15a) 

oz 

where  a is  the  coefficient  of  thermal  expansion  of  the  rod  material  (assumed  constant). 


Formally  integrating  this  expression,  one  obtains 
w{z,t)  = w(0,t)+af^T((,0d( 

J o 


iA-\5b) 


where  w(0,t)  is  the  thermal  expansion  of  the  rod  at  z = 0. 

We  are  interested  in  the  thermal  expansion  of  the  rod  relative  to  the  plane  z = 0.  Hence,  w(0,t)  = 
0.  Then,  to  compute  the  thermal  expansion  one  must  integrate  the  expression  for  the  temperature 
distribution  with  respect  to  the  distance  along  the  rod.  However,  it  is  first  necessary  to  define  a 
constant  reference  temperature  for  which  one  defines  the  state  of  zero  thermal  expansion.  That 
is,  the  temperature  change  of  the  rod  is  measured  relative  to  a defined  reference  temperature. 

Denoting  the  reference  temperature  as  Tcai  and  using  the  temperature  scale  given  by  Eq(A-6),  the 
expression  for  the  temperature  is: 

T(z,t)  = T„-T^,  + T*(z,t)  (A-16a) 

and  the  expression  for  the  thermal  expansion  is 

w(z,0  = cl{T^-TJz  + arr*(Ti,0^r|.  {A-\6b) 


The  explicit  expression  for  T*(z,t)  depends  upon  the  solution  of  the  temperature  problem  defined 
by  Eq(A-lO).  For  a constant  temperature  input,  the  expression  for  T*(z,t)  is  given  by  Eq(A-l  1). 
For  a constant  heat  flux  input,  the  expression  for  T*(z,t)  is  given  by  Eq(A-14).  Although  these 
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expressions  for  T*(z,t)  look  complicated,  they  can  be  integrated  to  give  an  explicit  functional 
expression  for  the  thermal  expansion.  These  expressions  are  presented  below.  First,  however,  a 
brief  discussion  of  the  above  expression  for  thermal  expansion,  Eq(A-16b),  is  worthwhile. 

First,  consider  that  the  rod  is  manufactured  with  a scale  along  its  length.  This  scale  may  be  a 
ruling  as  for  a common  ruler  or  it  may  be  a lead  screw.  Assume  that  the  rod  is  manufactured  and 
the  ruling  calibrated  at  a constant  temperature,  T^ai . (This  is  the  reason  for  the  nomenclature  T^aj 
for  the  reference  temperature.)  Now,  assume  that  the  rod  is  placed  in  an  environment  at  a 
temperature  T„  and  the  rod  experiences  a temperature  change  until  the  rod  is  at  the  steady 
temperature  T„.  This  temperature  change,  T„  - T^ai,  results  in  a thermal  expansion  of  the  rod 
which  is  given  by  the  first  term  in  Eq(A-16b).  This  thermal  expansion  is  relative  to  the  plane 
z = 0 of  the  rod.  The  thermal  expansion  results  in  a scale  error  that  is  directly  proportional  to  the 
temperature  change  and  the  calibrated  distance,  z,  from  the  end  of  the  rod.  At  the  uniform 
temperature,  T„,  the  measured  distance  using  the  manufactured  ruling  is: 

= z + w(z,0)  @ T.  (A- 17a) 

After  the  initiation  of  the  heating  of  the  rod  at  z = 0,  the  measured  distance  using  the 
manufactured  ruling  is: 

Zm=asured  = Z + w(Z,t).  (A- 17b) 

Now,  expressions  for  the  thermal  expansion  of  the  rod  are  obtained. 

A. 5. 2 Thermal  Expansion  for  Constant  Temperature  Input 

For  a constant  temperature  input,  the  temperature  response  T*(z,t)  is  given  by  Eq(A-l  1). 
Substituting  this  expression  into  Eq(A-16b)  and  integrating,  one  obtains  the  result 

w{z,t)  = a(r^-r^^^)z+av/K7X(r^-rj{sinh(zv^)+ery(yX/) 

-Mle  v/ia +vlf) 

+ l/2e  ■">/^ery(z/20a-v/Ir)}  {A  -18) 


where  T^  is  the  initial  temperature  input  to  the  rod  at  z = 0 (see  Eq(A-9)). 

The  limiting  values  for  the  expression  for  thermal  expansion  are: 

limit  w(z,t)  = 0 (A- 19a) 

z-0" 


limit  w(z,t)  = a(T„  - T^a,)  Zi^g,  + aV (k/X)  (To-T„)erf/ (Xt) 


(A- 19b) 


Note:  Zjarge  is  a large  value  of  z and  not  infinite 


limit  w(z,t)  = a(T„  - T^ai)  z (A- 1 9c) 

t-0" 

limit  w(z,t)  = a(T.  - T^,)  z + a/ (k/A)  (T,-T„)[  1 - (A- 1 9d) 

t-  oo 

The  first  limit  emphasizes  that  the  thermal  expansion  is  zero  at  the  z = 0 cross-section  of  the  rod. 
The  second  limit  gives  the  time  dependence  of  the  thermal  expansion  for  large  z.  The  third  limit 
is  the  initial  thermal  expansion  due  to  the  temperature  difference  between  the  calibration 
temperature  and  the  ambient  temperature.  The  last  limit  is  the  steady  state  (i.e.,  time 
independent)  thermal  expansion. 

A.5.3  Thermal  Expansion  for  Constant  Heat  Flux  Input 

For  a constant  heat  flux  input,  the  temperature  response  T*(z,t)  is  given  by  Eq(A-14). 
Substituting  this  expression  into  Eq(A-16b)  and  integrating,  one  obtains  the  result 

w{z,t)  = a(r^ +a(K/A)(^yA:){  1 -coshfryX/K) 

-e  ~^erf{zl2yfKt) 

+ 1 /2e  ^^erf{zl2sfKt+ylh) 

+ \/2e  ~^^er/{zl2y/Kt-yfXt)}  (A  -20) 

where  q^  is  the  initial  heat  flux  input  to  the  rod  at  z = 0 (see  Eq(A-12)). 

The  limiting  values  for  the  expression  for  thermal  expansion  are: 


limit  w(z,t)  = 0 
z-0" 

(A-21a) 

limit  w(z,t)  = a(T.  - T,„)  Z|„g,  + a/ (k/X)  (qo/k)[l  - e'^'] 

Z-  oo 

Note:  Ziarge  is  a large  value  of  z and  not  infinite 

(A-21b) 

limit  w(z,t)  = a(T„  - T^^i)  z 
t-0^ 

(A21c) 

limit  w(z,t)  = a(T„  - T^ai)  z + a/ (k/A)  (qo/k)[l  - 

t-  oo 

(A-21d) 
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The  first  limit  emphasizes  that  the  thermal  expansion  is  zero  at  the  z = 0 cross-section  of  the  rod. 
The  second  limit  gives  the  time  dependence  of  the  thermal  expansion  for  large  z.  The  third  limit 
is  the  initial  thermal  expansion  due  to  the  temperature  difference  between  the  calibration 
temperature  and  the  ambient  temperature.  The  last  limit  is  the  steady  state  (i.e.  time 
independent)  thermal  expansion.  The  reader  should  note  that  the  steady  state  thermal  expansion 
for  constant  heat  flux,  Eq(A-21d),  is  functionally  identical  to  the  steady  state  thermal  expansion 
for  constant  temperature  input,  Eq(A-19d). 
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APPENDIX  B - SURFACE  COOLING  DETAILS 


As  described  in  Appendix  A,  surface  cooling  along  the  length  of  the  rod  is  described  by  the  term 
A(T-TJ  in  Eq  (A-5b).  The  parameter  X defines  the  cooling  rate  and  is  expressed  in  terms  of  the 
physical  parameters  as 

X = hp/pcS  (B-1) 

where  h is  the  convection  heat  transfer  coefficient  between  the  lateral  surface  of  the  rod 
and  the  ambient  environment 
p is  the  perimeter  of  the  rod  cross-section 
S is  the  area  of  the  rod  cross-section 
p is  the  mass  density  of  the  rod  material 

c is  the  specific  heat  (at  constant  pressure)  of  the  rod  material. 

The  discussion  in  this  appendix  focuses  upon  the  estimation  of  the  heat  transfer  coefficient,  h. 
This  estimation  is  rather  “cookbook”  for  persons  familiar  with  heat  transfer  problems.  An 
outline  of  the  procedure  is  given  so  that  the  reader  may  appreciate  the  detail  considerations 
without  the  pain  of  a refresher  course  in  heat  transfer.  For  more  detail,  the  interested  reader  may 
consult  reference  [B-1]. 

B.l  Estimation  of  the  Convection  Heat  Transfer  Coefficient,  h 

For  practical  problems,  we  consider  the  rod  to  be  placed  in  air  at  a represenstative  room 
temperature  of  23  °C  (296.2  °K)  and  atmospheric  pressure.  The  properties  of  air  are  extensively 
tabulated  [B-1,  Table  A-5,  p.661]. 

First,  it  is  necessary  to  calculate  the  “film”  temperature  of  the  rod  surface-air  interface.  The  film 
temperature  is  defined  as 

Tf(z,t)  = (T(z,t)  + T.)/2,°K  (B-2) 

using  the  notation  of  Appendix  A on  the  right  hand  side  and  a subscript  f to  denote  “film”. 

The  film  temperature,  Tf,  is  used  to  look  up  the  properties  of  air  that  are  needed  to  estimate  the 
heat  transfer  coefficient,  h.  From  equation  (B-2)  we  see  that  the  film  temperature  is  a function  of 
both  position  along  the  rod  and  time.  Strictly  speaking,  the  heat  transfer  coefficient  also  is  a 
function  of  position  and  time.  Such  considerations  lead  to  non-linear  problems  in  heat 
conduction.  In  Appendix  A,  it  was  assumed  that  the  parameters  h and  X were  constants  so  that 
the  resulting  differential  equation  (Eq(A-7))  is  linear.  For  the  present  development  we  make  the 
simplifying  assumption  that  Tf=  T^  which  is  the  initial  temperature  condition  for  the 
developments  in  Appendix  A. 
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The  cooling  law  expressed  by  Eq(A-3)  corresponds  to  forced  convection  cooling  in  the  heat 
transfer  world.  The  heat  transfer  coefficient  is  determined  empirically  and  depends  upon  the 
cross-section  shape  and  the  orientation  of  the  transverse  air  velocity  relative  to  the  cross-section. 

Since  ball  screws  are  of  circular  cross-section,  the  following  discussion  is  limited  to  this  shape 
(i.e.,  our  rod  is  a long  circular  cylinder). 

Knowing  the  film  temperature  (Tf  = TJ  and  the  cross-section  shape,  one  is  then  in  a position  to 
compute  the  heat  transfer  coefficient.  For  a circular  cross-section  cooled  by  a transverse  low 
velocity  air  flow  the  Nusselt  number  is  empirically  related  to  the  Reynolds  number  and  the 
Prandtl  number  where  the  properties  of  air  are  evaluated  at  the  film  temperature  (see  Reference 
B-1,  Chapter  6).  The  Nusselt,  Reynolds  and  Prandtl  numbers  are  dimensionless  numbers 
comprising  groupings  of  the  physical  properties  of  the  fluid  surrounding  the  rod  and  a length 
parameter  taken  as  the  rod  diameter. 

For  the  present  problem  (or  notation)  we  have: 

Nusselt  number:  Nu  = hd/kf 

Reynolds  number:  Re  = du^Uf 

Prandtl  number:  Pr  = (pCpU/k)f 

where  d is  the  diameter  of  the  circular  cross-section  rod 
k is  the  thermal  conductivity  of  air 
u„  is  the  transverse  velocity  of  the  air 
p is  the  density  of  air 

Cp  is  the  specific  heat  of  air  at  constant  pressure 
V is  the  kinematic  viscosity  of  air 

subscript  f denotes  the  quantity  evaluated  at  atmospheric  pressure  and  Tf=T„ 

First,  one  must  calculate  the  Reynolds  number  based  on  the  rod  diameter.  For  the  present 
discussion,  we  consider  a range  of  rod  diameters,  cross-flow  air  velocity,  and  air  temperature. 
These  ranges  are  given  in  Table  B-1 . For  these  ranges,  the  range  of  Reynolds  number  is: 

190  < Re  < 13300 

and  the  range  of  Pandtl  number  is: 

0.689  < Pr  < 0.722. 

Knowing  the  Reynolds  number  and  the  Prandtl  number,  one  uses  an  empirical  relationship 
between  Re  and  Pr  to  estimate  the  Nusselt  number  for  the  rod-air  state.  Several  such 
relationships  are  available  and  are  usually  restricted  to  a specific  range  of  Reynolds  number.  For 
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example,  see  reference  [B-1]  Equation  (6-18).  For  the  range  of  parameters  described  above,  the 
range  of  Nusselt  number  is: 

8<Nu<71. 

The  importance  of  the  Nusselt  number  is  that  the  heat  transfer  coefficient,  h,  and  the  cooling 
coefficient,  A,  are  directly  proportional  to  the  Nusselt  number. 

Knowing  the  Nusselt  number,  one  immediately  calculates  the  heat  transfer  coefficient,  h.  The 
range  is: 

1.8  <h< 31.6,  W/m^X. 

To  calculate  the  cooling  coefficient,  A,  one  must  specify  the  rod  material  in  addition  to  the 
parameters  discussed  above.  As  presented  in  Table  B-1,  for  steel  rods  the  range  of  cooling 
coefficients  is  [B-2]: 

1.2- 10'' <X  <7.2*10-",  1/s. 

For  air  at  T„=  300 °K  = 26.85  °C  the  range  of  h and  X are: 

1.8  <h< 31.3  W/m2°C 

1.2- 10-' <X  <7.2-10-",  1/s. 

Hence,  one  sees  that  the  range  of  heat  transfer  coefficient,  h,  and  cooling  coefficient,  X,  for  a 
circular  cross-section  rod  is  essentially  a function  of  the  rod  diameter,  material  properties,  and  air 
velocity  and  is  not  strongly  dependent  upon  the  ambient  air  temperature. 

As  mentioned  above,  the  parameters  h and  X are  directly  proportional  to  the  Nusselt  number. 

The  explicit  expressions  are: 

h = kfNu/d  (B-3a) 

X = 4kfNu/(d^pc)  (B-3b) 

where  kf  is  the  thermal  conductivity  of  air  at  Tf  = T„ 

d is  the  diameter  of  the  circular  cross-section  rod 
pc  is  the  heat  capacity  of  the  rod  material. 
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Table  B-1 : Range  of  Parameters 


circular  cross-section  steel  rod  with  forced  convection  cooling  by  air  at  atmospheric  pressure 


material 

parameter 

range 

units 

air 

temperature 

250°  < T„  < 400° 

°K 

air 

Prandtl  number 

0.689  < Pr  < 0.722 

dimensionless 

air 

velocity 

0.1  <u„<  1.0 

m/s 

rod 

diameter 

0.05  <d<  0.15 

m 

rod 

ratio:  perimeter  to  area 

27  < p/s  <80 

1/m 

rod 

thermal  conductivity  of  steels 

1 1<  k < 60 

W/m°C 

rod 

heat  capacity  of  steels 

3.5-10^  < pc  <4.2-10^ 

rod 

thermal  diffusivity  of  steels 

M0-‘<k<2-I0-^ 

m^/s 

air-rod 

Reynolds  number 

190  < Re  < 13,300 

dimensionless 

air-rod 

Nusselt  number 

8<Nu<71 

dimensionless 

air-rod 

heat  transfer  coefficient 

1.8<h<31.6 

W/m^°C 

air-rod 

cooling  coefficient 

1.2-10-^<X<7.  2-10^ 

1/s 
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